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Abstract 

3-flavor quark matter (strange quark matter; SQM) can be stable or meta- 
stable for a wide range of strong interaction parameters. If so, SQM can play an 
important role in cosmology, neutron stars, cosmic ray physics, and relativistic 
heavy-ion collisions. As an example of the intimate connections between astro- 
physics and heavy-ion collision physics, this Chapter gives an overview of the 
physical properties of SQM in bulk and of small-baryon number strangelets; 
discusses the possible formation, destruction, and implications of lumps of SQM 
(quark nuggets) in the early Universe; and describes the structure and signa- 
ture of strange stars, as well as the formation and detection of strangelets in 
cosmic rays. It is concluded, that astrophysical and laboratory searches are 
complementary in many respects, and that both should be pursued to test the 
intriguing possibility of a strange ground state for hadronic matter, and (more 
generally) to improve our knowledge of the strong interactions. 



1 Introduction 



Hadronic matter is expected to undergo a transition to quark-gluon plasma under 
conditions of high temperature and/or baryon chemical potential. These conditions 
may be achieved for a brief moment in ultrarelativistic heavy-ion collisions, but they 
are also likely to appear in Nature. A very high density (and comparatively low 
temperature) environment exists in the interior of neutron stars, which may actually 
contain significant amounts of quark matter in the interior. High temperatures (but 
rather low baryon chemical potential) were realized in the first 10~ 4 seconds after 
the Big Bang, and here a hot quark-gluon plasma state must have existed until the 
temperature dropped to 100-200 MeV due to the adiabatic expansion of the Universe. 

This Chapter will outline some of the possible ways in which astrophysics may 
teach us about the existence and properties of quark-gluon plasmas. The advantage 



1 



relative to laboratory searches is, that truly bulk systems can be studied, and that the 
timescales involved are much longer than those relevant to collisions. Disadvantages 
are that astrophysicists (with possible exceptions if strange quark matter is absolutely 
stable) can only observe indirect consequences of the plasma state for example in the 
properties of pulsars or in the distribution of light nuclei produced a few minutes 
after the Big Bang. It will be shown, however, that astrophysics arguments in many 
cases can be used to constrain parameters significantly relative to direct experimental 
approaches because of the large volumes and timescales involved. 

The implications of quark-gluon plasmas in astrophysics and cosmology are many- 
fold, and I shall focus on aspects related to the idea of (meta) stability of strange quark 
matter through discussions partly biased by my own research interests. 

Lumps of up, down, and strange quarks (strange quark matter, SQM), with masses 
ranging from small nuclei to neutron stars, rather than 56 Fe, could be the ground state 
of hadronic matter even at zero temperature and pressure. This possibility, first noted 
by Bodmer in 1971 [Q, has attracted much attention since Witten resurrected the 
idea in 1984 ||. The existence of stable or metastable SQM would have numerous 
consequences for physics and astrophysics, and testing some of these consequences 
should ultimately tell us whether SQM really exists. 

First it was believed that SQM might give a natural explanation of the cosmo- 
logical dark matter problem. While not ruled out, this idea is now less popular, but 
strange quark matter may still be important in astrophysical settings, such as strange 
stars. Numerous investigations have searched for deposits of SQM on the Earth and 
in meteorites, so far unsuccessfully, and recently relativistic heavy-ion collision exper- 
iments have been performed and/or proposed to test the idea. Cosmic ray searches 
have come up with a few potential candidates for small SQM-lumps (strangelets), but 
at present no compelling evidence for stable SQM has been presented. This, however, 
does not rule it out. Most searches for SQM are sensitive to strangelets with very 
low baryon number, A, and as discussed later, finite size effects have a significant 
destabilizing effect on such objects, even if SQM is stable in bulk. 

There is a significant range of strong interaction parameters for which SQM in 
bulk is stable. But even if it is not, many of the (astro)physical implications are more 
or less unchanged in the case of metastable SQM. In neutron stars, for instance, the 
high pressure brings SQM closer to stability relative to hadronic matter, and it is 
quite likely, that neutron stars contain cores of strange quark matter, even if SQM is 
unstable at zero pressure. In relativistic heavy-ion collision experiments, strangelets 
need "only" survive for 10 -8 seconds to be of interest. In fact, (meta) stable strangelets 
may be one of the "cleanest" signatures for formation of a quark-gluon plasma in such 
collisions. 

The present review tries to give an account of the status of strange quark matter 
physics and astrophysics, as of early 1998, but of course not all aspects are covered 
in equal detail. In particular, nothing is said about the heroic experimental efforts 
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to produce strangelets in heavy-ion collisions. A collection of papers describing all 
aspects of SQM and a list of references to the field through mid-1991 can be found in 
H • An earlier review was given in [[| . Recent reviews include |5], ||, 0, [|, £| [10], [11], [12 



and the thorough reader will notice, that some parts of the present Chapter borrows 
from my own papers among these since the physics discussed therein remains more 



or less unchanged. Refs. [11, 12| also discuss the related issue of lumps of metastable 



strange hadronic matter, which will not be dealt with here. 

Section |2| discusses the physics of SQM, starting out with simple estimates of why 
3-flavor quark-matter is likely to be more bound than the 2-flavor alternative, proceed- 
ing with more detailed descriptions of SQM in bulk. Smaller systems (strangelets), 
for which finite-size effects are crucial, are described in Section ||. Most of the results 
are based on the MIT bag-model, but it is worth stressing from the outset that this 
should only be viewed as a crude approximation to reality, ultimately to be surpassed 
by direct QCD-calculations. 

Section [| deals with the possible production of lumps of SQM (often called quark 
nuggets) in the cosmological quark-hadron phase transition, and the struggle of quark 
nuggets to survive evaporation and boiling in a hostile environment. It turns out, 
that only large nuggets are likely to survive, but the physics involved in the destruc- 
tion process is illuminating as it resembles the (time-reversed) physics involved in 
strangelet production in heavy-ion collisions. Implications of surviving quark nuggets 
for Big Bang nucleosynthesis and the dark matter problem are also discussed. 

Perhaps the most likely place to discover SQM (even if it is not absolutely stable) 
is in neutron stars. These could be "hybrid", "strange", or even "mixed" (the first 
term conventionally used for neutron stars with quark cores; the second for "true" 
quark stars in case of SQM stability, and the latter for objects with mixed phases of 
quark matter and nuclear matter). Section |5] describes these stars, their implications 
for our understanding of pulsars, and the possible connection to the energetic gamma- 
ray bursters. 

Strangelets surviving from the early Universe or released from strange stars in 
binary systems have been searched for in cosmic ray detectors and in meteorites and 
mineral deposits. So far there are only a few potential candidates, but more sensitive 
experiments will soon be carried out. Section ^] discusses some of the limits obtained. 
It also presents an astrophysical argument which either improves the Earth-based 
flux-limits by many orders of magnitude (almost excluding absolutely stable SQM), 
or predicts that all neutron stars are strange stars, if SQM is stable (the prediction to 
choose depends on whether any pulsars can be proven to be ordinary neutron stars). 

Conclusions and a brief outlook are provided in Section [j]. 



2 Physics of SQM in bulk 



3 



2.1 Does strange matter conflict with experience? 



At first sight, the possibility that quark matter could be absolutely stable seems 
to contradict daily life experiences (and experiments) showing that nuclei consist of 
neutrons and protons, rather than a soup of quarks. If a lower energy state exists, 
then why are we here? Why have we not decayed into strange quark matter? 

The answer to this obvious question is, that (meta) stability of strange quark 
matter requires a significant fraction of strange quarks to be present. Conversion 
of an iron nucleus into an A = 56 strangelet thus demands a very high order weak 
interaction to change dozens of u- and <i-quarks into s-quarks at the same time. 
Such a process has negligible probability of happening. For lower A the conversion 
requires a lower order weak interaction, but as demonstrated later, finite-size effects 
destabilize small strangelets so that they become unstable or only weakly metastable 
even if strange quark matter is stable in bulk. 

Therefore (meta) stability of strange quark matter does not conflict with the exis- 
tence of ordinary nuclei. On the other hand, the existence of ordinary nuclei shows, 
that quark matter composed of u- and <i-quarks alone is unstable, a fact that will be 
used later on to place constraints on model parameters. 

Another constraint from our mere existence can be placed on the electrical 
charge of strangelets. If energy is gained by converting ordinary matter into strange 
quark matter, strangelets with negative quark charge, even if globally neutral due to 
a cloud of positrons, would have devastating consequences, eating up the nuclei they 
would encounter. Even a small stable component in the cosmos would be intolerable 
(but they could still appear as metastable products in heavy-ion collisions, like the 
recent charge —1, mass 7.4 GeV event in NA52 at CERN fl~3||). A positive charge on 
the quark surface (neutralized by surrounding electrons) is less problematic, because 
ordinary nuclei will be electrostatically repelled. The barrier has to be of a certain 
height, though, in order not to impact stellar evolution (see below). Note that neu- 
trons are easily absorbed. As demonstrated later, this has important consequences 
for quark star formation and can be used to constrain strange matter properties using 
several astrophysical lines of reasoning. It may even lead to practical applications in 



energy production, etc. [14 



2.2 Simple arguments for (meta) stability 

As argued above, quark matter composed of u and d-quarks is expected to be unstable 
(except from 3-quark baryons). Introducing a third flavor makes it possible to reduce 
the energy relative to a two-flavor system, because an extra Fermi-well is available. 
The introduction of an extra fermion-flavor makes it possible to increase the spatial 
concentration of quarks, thereby reducing the total energy. A penalty is paid because 
the mass of the s-quark is high compared to that of u and d, so stability is most likely 
for low s-quark mass. 
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To make the argument slightly more quantitative, consider non-interacting, mass- 
less quarksQ inside a confining bag at temperature T = 0, without external pressure. 
For a massless quark-flavor, i, the Fermi momentum, ppi, equals the chemical po- 
tential, fii (throughout the chapter, unless otherwise noted, h — c — ks — 1; for 
an introduction to Fermi-gas thermodynamics, see for instance Ref. |15[). Thus the 
number densities are rii = fi^/n 2 , the energy densities e« = 3/z 4 /(47r 2 ), and the pres- 
sures Pi = /z 4 /(47r 2 ). The sum of the quark pressures is balanced by the confining bag 
pressure, B; J2i Pi — B; the total energy density is e = J2i e % + B = 3 J2i Pi + B = 4B, 
and the density of baryon number is n B = J2i n i/3. Notice that the sum of the con- 
stituents pressures, as well as the total energy density are given solely in terms of the 
bag constant, B. 

For a gas of u and d-quarks charge neutrality requires rid = 2n u , or ii 2 = 
H u = 2~ 1 / 3 n c [- The corresponding two-flavor quark pressure is P 2 = P u + Pd = 
(1 + 2 4//3 )/i2/(47r 2 ) = B, the total energy density e 2 = 3P 2 + B = AB, and the baryon 
number density n B2 = (n u + n d )/3 = /u|/vr 2 , giving an energy per baryon of 

e 2 /n B2 = (1 + 2 4 / 3 ) 3 / 4 (4tt 2 ) 1 / 4 j B 1 / 4 = 6.441E 1 / 4 w 934MeV J B 1 1 4 / 5 4 , (1) 

where B\^ = 5 1 / 4 /145MeV; 145MeV being the lowest possible choice for reasons 
discussed below. 

A three-flavor quark gas is electrically neutral for n u = rid = n s , i. e. /i 3 = fi u = 
— F° r fixed bag constant the three-quark gas should exert the same pressure as 
the two-quark gas (leaving also the energy density, e 3 = 3P 3 + B — 4B, unchanged). 
That happens when /i 3 = [(1 + 2 4 / 3 )/3] 1//4 /i2, giving a baryon number density of 
n B3 = /^lY 71 " 2 = [(1 + 2 4//3 )/3] 3 / 4 ne2- The energy per baryon is then 

e 3 /n m = 3/i3 = 3 3/4 (4vr 2 ) 1 / 4 J B 1 / 4 = 5.714P 1 / 4 « 829MeV BlU; (2) 

lower than in the two-quark case by a factor n B 2/n B3 = (3/(1 + 2 4 / 3 )) 3 / 4 w 0.89 . 

The possible presence of electrons was neglected in the calculations above. For 
two-flavor quark matter, including electrons in chemical equilibrium via u + e~ «-> 
d + v e , so that n u + \x e = fid, gives more cumbersome equations, but only changes 
C2,Mb2 to 6.445-B 1 / 4 , since \i e turns out to be rather small. Three-flavor quark matter 
does not contain electrons for non-interacting, massless quarks. 

One may therefore gain of order 100 MeV per baryon by introducing an extra 
flavor. At fixed confining bag pressure the extra Fermi- well allows one to pack the 
baryon number denser into the system, thereby gaining in binding energy. 

The energy per baryon in a free gas of neutrons is the neutron mass, m n = 
939.6MeV; in a gas of 56 Fe it is 930 MeV. Naively, stability of w<i-quark matter rel- 
ative to neutrons thus corresponds to e 2 /nB2 < Tn n , or B 1 / 4 < 145.9MeV (5 1 / 4 < 

1 Since current quark masses rather than constituent quark masses enter in the MIT bag model 
used to describe SQM, this is a very good approximation for u and d-quarks with 5 MeVw m u < 
m d w 10 MeV < 300 MeV « fj, u , fi d . 
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144.4MeV for stability relative to iron). The argument can be turned around: Since 
one observes neutrons and 56 Fe in Nature, rather than ud-quark matter, it is con- 
cluded that B 1 / 4 must be larger than the numbers just quoted. More detailed calcu- 
lations including finite-size effects and Coulomb-forces do not change these numbers 
much, so we shall assume for the present purpose that 5 1/4 = 145MeV is an experi- 
mental lower limit for a s = 0. (Here a s denotes the strong "fine-structure" constant; 
a s = corresponding to non-interacting quarks except for the confinement given by 
B)- 

Bulk strange quark matter is absolutely stable relative to a gas of iron for B 1 ! 4 < 
162.8MeV, metastable relative to a neutron gas for B 1 / 4 < 164.4MeV, and relative 
to a gas of A-particles (the ultimate production limit in heavy-ion collisions) for 
B 1 / 4 < 195.2MeV. These numbers are upper limits. As demonstrated below, a finite 
s-quark mass as well as a non-zero strong coupling constant decreases the limit on 
B 1 ' 4 . 

The presence of ordinary nuclei in Nature cannot be used to turn the values of 
B 1 ! 4 just quoted for SQM into lower limits. Conversion of a nucleus into a lump 
of SQM requires simultaneous transformation of roughly A u- and <i-quarks into s- 
quarks. The probability for this to happen involves a weak interaction coupling to 
the power A, i. e. it does not happen. This leads to the conclusion, that even if 
SQM is the lowest energy state for hadronic matter in bulk, its formation requires 
a strangeness-rich environment or formation via a "normal" quark-gluon plasma in 
relativistic heavy-ion collisions, the early Universe, or a neutron star interior. All of 
these possibilities will be explored in the following. 



2.3 SQM in bulk at T = 

The estimates above assumed m s = a s = 0. Non-zero a s was found by Farhi and Jaffe 



lq] to correspond effectively to a reduction in B. In the interest of simplicity I will 
therefore set a s = in most of the following. The energy "penalty" paid by having 
to form s-quarks at a finite mass of 50-300 MeV calls for more detailed calculations, 
however. Such calculations are usually performed within the MIT bag model |L7|, [18| . 

Strange quark matter contains degenerate Fermi gases of u, d, and s quarks, and 
e~ or e + . Chemical equilibrium is maintained by weak interactions, 

d *—> u + e~ 4~ v e (3) 
s «-> u + e~ + v e (4) 
u + s <-> d + u, (5) 

where the first two reactions should be understood to include also the various per- 
mutations of the involved particles. 

Neutrinos generally escape the system, so we shall ascribe to them no chemical 



6 



potential. Thus the chemical potentials in equilibrium are given by 

/id = Us = Vu + He- (6) 

Knowing the chemical potentials one can calculate the thermodynamic potentials. 

(7) 
(8) 
(9) 
(10) 
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defining A = m s j\i s . 

Number densities are given by 

Uiy = -dfliy/dfif, (11) 

i. e. n e y = /ie/3vr 2 , n u y = /i„/vr 2 , n d y = nl/ir 2 , n s y = /i^(l - A 2 ) 3/2 /7r 2 . The total 
pressure is 

P = Y f P i -B = - y £n itV -B = 0, (12) 

i i 

and charge neutrality requires 

2 11 

2 n u,v ~ ^n d y - -n s y - n e y = 0. (13) 
The total energy density is 

e = J2(tt i ,v + n i y{i i ) + B, (14) 

i 

and the density of baryon number 

n B = ^(n u y + n dy + n s y). (15) 

Combining Eqs. (§) and (13) leaves only one independent chemical potential, 
which can be determined from the pressure balance, Eq. (0). Thus all parameters 
can be calculated for a given choice of m s and B. Results of such calculations are 
shown in Figure [I]. Similar calculations were originally done by Farhi and Jaffe [T(| . 

The calculations above assumed zero temperature and external pressure. Finite 
temperature and external pressure can be relevant in connection with cosmology 
(Section |j) and strange stars (Section |]) respectively, and also for strangelet creation 
in collision experiments. The relevant extensions of the formulae above will be given 
in Section 
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Figure 1: Energy per baryon for bulk strange quark matter as a function of bag 
constant and strange quark mass. 

3 Strangelets 

So far the treatment of SQM has focused on the bulk properties. This approximation 
is generally valid for large baryon numbers. For A <C 10 7 the quark part of SQM is 
smaller than the Compton wavelength of electrons, so electrons no longer ensure local 
charge neutrality. Therefore Coulomb energy has to be taken into account, though 
the fortuitous cancellation of q u +Qd+Qs = § — | — § = ^ means that Coulomb energy is 
much less important for strangelets than for nuclei. For even smaller baryon numbers 
(in practice A < 10 3 ) other finite size effects such as surface tension and curvature 
have to be taken into account. 

Several strangelet searches with relativistic heavy-ion collisions as well as cosmic 
ray searches have been carried out, and others are planned for the future. Most 
of these searches are sensitive only to low A-values, so it is important to know the 
properties of small lumps of strange quark matter (strangelets). 

In the following I will describe the physical properties of strangelets in the lan- 
guage of the MIT-bag model (only limited work has been performed using other 
models — qualitatively confirming the MIT-bag results, though quantitative details 
can differ). First, I will discuss results obtained from direct solution of the Dirac 
equation with MIT-bag boundary conditions; such mode filling calculations corre- 
spond to a nuclear shell model. Then I will show how the mean behavior of the shell 
model results can be understood physically in terms of a liquid drop model calculation 
based on a smoothed density of states, and how approximations to the liquid drop 
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results give simple formulae for strangelet masses etc. Finally I discuss the changes 
introduced if strangelets are at finite rather than zero temperature. 



3.1 Shell model 

Mode-filling for large numbers of quarks in a spherical MIT-bag was performed for 
u<i-systems by Vasak, Greiner and Neise ||19|| , and for 2- and 3-flavor systems by 
Farhi and Jaffe |L6[], and Greiner et al. |2(| (see also pi]). Gilson and Jaffe |22| 



published an investigation of low-mass strangelets for 4 different combinations of s- 
quark mass and bag constant with particular emphasis on metastability against strong 
decays. Further parameter ranges were studied and compared to liquid drop model 
calculations by Madsen , and recently new shell-model studies were published by 



Schaffner-Bielich et al. [£4j]. All of these calculations were performed for a s = 0, 
which will also be assumed in the following. 

In the MIT bag model noninteracting quarks are confined to a spherical cavity of 
radius R. They satisfy the free Dirac equation inside the cavity and obey a boundary 
condition at the surface, which corresponds to no current flow across the surface. 
The bag itself has an energy of BV . In the simplest version the energy (mass) of the 
system is given by the sum of the bag energy and the energies of individual quarks, 

E= E E^K 2 + <i) V2 + S4 ^ 3 /3. (16) 

i=u,d,s ft 

Here k K i = x Ki i/R, where X i £1X6 61 genvalues of the equation 

f ^ K ' l) = «■ + m\R^ + mi R fK -^ XK ' l) - (17) 
f K are regular Bessel functions of order k , 

f ( r ) — I j*( X ) (1Q) 

U[X) \ y K (x) = (-l) K+1 i-«-i(x) k<0 [L * } 

For states with quantum numbers (j, /) k takes the values k = ±(j + |) for / = j ± |. 
For a given quark flavor each level has a degeneracy of N Kti = 3(2j + 1) (the factor 3 
from color degrees of freedom). For example, the IS1/2 ground-state (j = 1/2, / = 0, 
K — — 1) for a massless quark corresponds to solving the equation tanx = x/(l — x), 
giving x ~ 2.0428. The ground state has a degeneracy of 6 per flavor. 

For massless quarks (finding the equilibrium radius from dE/dR = 0) one gets 

E = 364.00MeV J B 1 1 4 / 5 4 x K ^' A (19) 

where the sum is to be taken over all 3A quark-levels, and the numbers x K ^ for 
massless quarks are tabulated in [|Hf|. 
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For massive quarks the level filling scheme is more cumbersome (see e.g. Refs. 
2^, ^3|). Fixing bag constant and quark- masses, for each baryon number one must 
fill up the lowest energy levels for a choice of radius; then vary the radius until a 
minimum energy is found (dE/dR = 0). Since levels cross, the order of levels is 
changing as a function of R. This is easily seen in the Figures, where one notices 
discontinuous changes in the position of shells. 
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Figure 2: Energy per baryon (in MeV) for strangelets with B 1 ^ = 145MeV and m s 
from 0-300 MeV in steps of 50 MeV (m s increases upward). The figure on the right 
shows an expanded view of the low-mass region to highlight the change of "magic 
numbers" with changing m s . 

One notices that the energy per baryon smoothly approaches the bulk limit for 
A — > oo, whereas the energy grows significantly for low A. For low s-quark mass 
shells are recognized for A = 6 (3 colors and 2 spin orientations per flavor), and 
less conspicuous ones for A = 18, 24, 42, 54, 60, 84, 102 etc. As m s increases it 
becomes more and more favorable to use u and d rather than s-quarks, and the 
"magic numbers" change; for instance the first closed shell is seen for A = 4 rather 
than 6. 

Equation (|16|) can be modified by inclusion of Coulomb energy and zero-point 
fluctuation energy. As already discussed the Coulomb energy is generally small. 
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The zero-point energy is normally included as a phenomenological term of the form 
—Z Q /R, where fits to light hadron spectra indicate the choice Z Q = 1.84. This was 
used, for instance, by Gilson and Jaffe ||22|| . Roughly half of this phenomenological 
term is due to center-of-mass motion, which can be included more explicitly by sub- 
stituting (I]x Kj j) 2 — Y,x 2 K i instead of (J2 x Kj i) 3 ^ in Eq. (]T9|) . The proper choice 



of a s and Z is a tricky question. As discussed by Farhi and Jaffe |16| the values are 
intimately coupled to B and m s , and it is not obvious that values deduced from bag 
model fits to ordinary hadrons are to be preferred. This uncertainty may have an 
important effect for A < 5-10, but the zero-point energy quickly becomes negligible 
for increasing A for reasons explained in Section |3.2.2| . It means, however, that it is 
difficult to match strangelet calculations to experimental data concerning ordinary 
hadrons or limits on the putative A = 2 iJ-dibaryon. 

3.2 Liquid drop model 

Mode-filling calculations are rather tedious but do of course give the "correct" results 
as far as the model can be trusted. But for many applications a global mass-formula 
analogous to the liquid drop model for nuclei is of great use and also gives further 
physical insight. 

A phenomenological approach to a strangelet mass-formula was undertaken by 
Crawford et al. f25|, E1J whereas Berger and Jaffe [27] made a detailed analysis within 



the MIT bag model. They included Coulomb corrections and surface tension effects 
stemming from the depletion in the surface density of states due to the mass of the 
strange quark. Both effects were treated as perturbations added to a bulk solution 
with the surface contribution derived from a multiple reflection expansion. Mad- 
sen [23, 28, 29] gave a self-consistent treatment including also the very important 



curvature energy. 

The following discussion closely follows [23|. All calculations are done for zero 
temperature and strong coupling constant, ot s . As argued by Farhi and Jaffe (16 
the latter assumption can be relaxed by a re-scaling of the bag constant. Also, I 
shall concentrate on systems small enough (A < 10 7 ) to justify neglect of electrons. 
Strangelets with A -C 10 7 are smaller than the electron Compton wavelength, and 
electrons are therefore mainly localized outside the quark phase. Thus strangelets do 
not obey a requirement of local charge neutrality, as was the case for SQM in bulk. 
This leads to a small Coulomb energy, which is rather negligible for the mass-formula 
(less than a few MeV per baryon), but which is decisive for the charge-to- mass ratio of 
the strangelet. A characteristic of strangelets, which is perhaps the best experimental 
signature, is that this ratio is very small compared to ordinary nuclei. Finally, I 
neglect charge screening, an issue of negligible importance for the mass formula, but 
of some importance for the charge-to-mass ratio for systems of radii above 5-10 fm 

(a > io 2 -io 3 ) pg. 
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In the ideal Fermi-gas approximation the energy of a system composed of quark 
flavors % is given by 

E = J2(^i + Nifn) + BV + E Coul . (20) 

i 

Here Qi, Ni and fii denote thermodynamic potentials, total number of quarks, and 
chemical potentials, respectively. B is the bag constant, V is the bag volume, and 
-E-Coui is the Coulomb energy. 



In the multiple reflection expansion framework of Balian and Bloch | 5T| , the ther- 
modynamical quantities can be derived from a density of states of the form 



dNi 
~dk 



.{£ + /.(=9« + /c(?i 0+ -.. 



(21) 



where area S = § dS (= AttR 2 for a sphere) and extrinsic curvature C = § (j^ + -^J dS 
(= 8ttR for a sphere). Curvature radii are denoted R\ and R 2 - For a spherical system 
Ri = R 2 = R- The functions fs and fc will be discussed below. 

In terms of volume-, surface-, and curvature-densities, n^y, n^s, and n^c, the 
number of quarks of flavor i is 



N 



dk 



dk = riiyV + Ui.sS + rii.cC. 



i,S>- 



(22) 



with Fermi momentum k Fi = (/i? — ml) 1 / 2 = //j(l — A 2 ) 1 / 2 ; Aj = mi/fa. 
The corresponding thermodynamic potentials are related by 



where cK^/cfyi; 



-Ni, and dQij/dfii = —Uij. The volume terms are given by 



n 



i,V 



4tt 2 



[1 - AD 1/2 (1 - l^) + ^ln 1 + (1 A A - )1/2 



rii,v 



^(1-A, 2 ) 3/2 . 



The surface contribution from massive quarks is derived from 



as [27| 
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(24) 
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(26) 
(27) 
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(1-Af) 1 



— I tan 

7T 



-1 



A,; 



-A.fl-A 2 ) 1 / 2 



(28) 



For massless quarks Qi t s = n^s = 0, whereas /c(0) = — 1/247T 2 gives 
l\ c = ^/8vr 2 ; n iiC = -/ii/4vr 2 . 

The curvature terms have never been derived for massive quarks, but as shown 
by Madsen [23|, the following Ansatz (found from analogies with the surface term 
and other known cases) works: 



*l7 




(29) 



This expression has the right limit for massless quarks (fc = —l/24n 2 ) and for 
infinite mass, which corresponds to the Dirichlet boundary conditions studied by 
Balian and Bloch [RTJ (fc = 1/127T 2 ). Furthermore, the expression gives perfect fits 
to mode-filling calculations (see the Figures and discussion below). From this Ansatz 
one derives the following thermodynamical potential and density: 



n 



i,C 



A? log 



a - \ 2 ) i/2 

A,; 



7T 



37rAj 



+ 7n 7^ 1" ttA, — — tan 



(1 " A 2 ) 1/2 



2A 



A, 



A; 



(30) 



riLC 



2 Xi Aj Aj 
With these prescriptions the differential of E(V, S, C, iVj) is given by 
dE = J2 (ttiydV + Q i>s dS + VL iyC dC + ^dJV<) + £dV + rfE Coul . 



(31) 



(32) 



Minimizing the total energy at fixed iVj by taking dE = for a sphere gives the 
pressure equilibrium constraint 



2 2 



^-^Coul 



(33) 



with 



Eqou\ 

dE C oui 



aZ v aZ 2 

Tor + ~2R' 

aZl aZ 2 



(34) 
(35) 



dV AOtcR 4 87T.R 4 ' 

where Zy = J2i QiTkyV is the volume part of the total charge, Z, whereas charge 
Z — Zy = J2i <li( n i,sS + n i,cC) is distributed on the surface. The quark charges are 
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q u = 2/3, qd = q s = —1/3. Eliminating B from Eq. ( p0|) then gives the energy for a 
spherical quark lump as 

E = £(iY>, + hl ijS S + \n ifi C) + \e Co vx. (36) 

■ 6 6 6 

The optimal composition for fixed baryon number, A, can be found by minimizing 
the energy with respect to Ni at fixed V, S, and C giving 



= dE = E 



dE, 



Coul 



dN,. 



dNi. (37) 



3.2.1 Massless quarks — Bulk limit 



For uncharged bulk quark matter Eq. ( p6|) reduces to the usual result for the energy 
per baryon 

e° = A-'Y.N^l (38) 

i 

where superscript denotes bulk values. The energy minimization, Eq. (|33|) , corre- 
sponds to 

* = -£^v = E^- (39) 

i i 

The last equality assumes massless quarks. In the bulk limit the baryon number 
density is given by 

«H?^- (40) 

and one may define a bulk radius per baryon as 

R° = (3/4t™° ) 1/3 . (41) 

For quark matter composed of massless u, d, and s-quarks, the Coulomb energy 
vanishes at equal number densities due to the fact that the sum of the quark charges 
is zero. Thus it is energetically most favorable to have equal chemical potentials 
for the three flavors. From the equations above one may derive the following bulk 
expressions for 3-flavor quark matter: 



'An 2 B' 



1/4 



1.905E 1 / 4 = 276.2MeV J B 1 1 £; (42) 



n° A = (^) 3 /tt 2 = 0.7005 3 / 4 (43) 
R° = (3/47r<) 1/3 = 0.699£~ 1/4 . (44) 
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And the energy per baryon is 



e ° = 3^° = 5.714B 1 / 4 , 



(45) 



in agreement with Eq. @. 

Following Berger and Jaffe []27[ one may to first order regard Coulomb, surface 
(and here correspondingly curvature) energies as perturbations on top of the bulk 
solution. In this approach one gets 



E 
A 



Q13/12 Dl/4 
e + A 2^ iZ i,C U - 6 + ^1/691/ 



./6 2 l/6^2/3 



829MeV + 351MeVA~ 2/3 



5 



1/4 
145- 



The corresponding result for 2-flavor quark matter (c.f. 
E 



is 



— = e° + A- 1 ^lc c ° ~ 934MeV + 291MeVyT 2/3 



B 



1/4 
145- 



(46) 



(47) 



3.2.2 Massive s-quarks — Bulk limit 

For m s > the energy minimization, Eq. fl39|), changes to 

b = -£n° v 



E 



+ 



, 47T 2 47T 2 

and the baryon number density is now given by 



1 - A 2 ) 1/2 (l - -A 2 ) + -A 4 In 1 + (1 7 A2)V2 
' v 2 ' 2 A 



1 



J2 %1 + %L(1-A 2 ) 3 / 2 



i=u,d 



7T 



7T- 



(48) 



(49) 



A bulk radius per baryon is still defined by Eq. QUI) . 

In bulk equilibrium the chemical potentials of the three quark flavors are equal, 
/i° = n° d = /j,® = fi° = e°/3. Neglecting Coulomb energy one may approximate the 
energy per baryon of small strangelets as a sum of bulk, surface and curvature terms, 
using the chemical potential calculated in bulk: 



E 
A 



(50) 



where S° = 4tt( j R°) 2 A 2 / 3 and C° = 8vr( J R°)A 1 / 3 . Examples for B 1 / 4 = 145MeV are 
(with s-quark mass in MeV given in parenthesis) 
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e(0) = 829MeV + 0MeVA- 1/3 + 351MeVA- 2/3 (51) 

e(50) = 835MeV + 61MeVA- 1/3 + 277MeVA~ 2/3 (52) 

e(150) = 874MeV + 77MeVA" 1/3 + 232MeVA" 2/3 (53) 

e(200) = 896MeV + 53MeVA" 1/3 + 242MeVA" 2/3 (54) 

e(250) = 911MeV + 22MeVA" 1/3 + 266MeVA~ 2/3 (55) 

e(300) = 917MeV + 0.3MeVA~ 1/3 + 295MeWT 2/3 (56) 

e(350) = 917MeV + 0MeVA- 1/3 + 296MeWT 2/3 (57) 



1050 




850 1 i i i i I i i i i I i i i 
50 100 

A 

Figure 3: Shell-model and liquid drop model results compared for B 1 ^ = 145MeV 
with massless u and d quarks, and with m s in the range 50-300 MeV in steps of 
50 MeV. For each value of m s the upper smooth curve is the full liquid drop model 
result, whereas the lower smooth curve is the bulk approximation. 

The bulk approximations above generally undershoot the correct solution with 
properly smoothed density of states by 2MeV for A > 100, 5MeV for A w 50, lOMeV 
for A w 10 and 20MeV for A w 5 (Figure §). This is because the actual chemical 
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potentials of the quarks increase when A decreases, whereas the bulk approxima- 
tions use constant \i. For massless s-quarks the expression for e(0) scales simply as 
B l /\ The same scaling applies for m s > e°/3, where no s-quarks are present; in 
the example above the scaling can be applied to e(350). For intermediate s-quark 
masses both e and m s should be multiplied by to scale the results. For in- 
stance, if B 1 ' 4 = 165MeVone finds e(150) = 985MeV + 93MeVA- 1 / 3 + 265MeVA- 2 / 3 ; 
e(250) = 1027MeV + 46M6VA' 1 / 3 + 284MeVA" 2 / 3 . Coulomb effects were not in- 
cluded above. Their inclusion would have no influence for m s — > 0, but would 
change the results by a few MeV for large m s . In particular charge neutral ud- 
quark matter has e = [934MeV + 291MeVA- 2 / 3 ] B\{1 (Eq. ©) rather than the 

917MeV + 296MeVA- 2 / 3 ] S}£ found above (Eq. (0)). 

In connection with the shell-model calculations I described the effects of a zero- 
point energy of the form —Z /R, and claimed that it was important only for A < 10. 
This can be understood in the bulk approximation of constant fi, because the zero- 
point term per baryon is proportional to A™ 4 / 3 compared to A -1 / 3 and A~ 2 I 3 for 
surface and curvature energies. The full term to be added to the bulk approximation 
expressions for a given e° is: 

e zero = -Z (4/243tt) 1 / 3 [ 2 + [1 - (3m s /e )] 3 / 2 ] 1/3 e°A~ 4 / 3 , (58) 

typically of order -200Z MeVA~ 4 / 3 . 



3.3 Shell model versus liquid drop model 

Self-consistent solutions can be obtained from Eq. (|36|) . These solutions are compared 
to the shell-model calculations and the bulk approximations in the Figures. The fits 
are very good, showing that inclusion of surface tension and curvature energy via the 
multiple reflection expansion explains the overall behavior of the results. 

3-flavor quark matter is energetically favored in bulk, and could be absolutely 
stable relative to 56 Fe for 144MeV < B 1 ^ < 163MeV. The lower limit corresponds 
to experimentally excluded stability of ud quark matter, whereas the upper limit 
corresponds to a bulk energy per baryon of w<is-matter of 930 MeV for m s = 0. 

Finite-size systems are strongly destabilized by the curvature energy, with a mag- 
nitude of about 300MeVyl _2 / 3 i?i45 for 3 quark flavors. This may pose problems for the 
experimental attempts of producing strange quark matter, since these experiments 
so far can only hope to create quark lumps with baryon number A < 20-30, and 
observe lifetimes exceeding 10~ 8 seconds. Further destabilization occurs for finite- 
mass s-quarks, where the surface tension (exactly zero for massless quarks) adds up 
to dOMeVA- 1 / 3 to the energy. 

Writing E/A = e° + c suli A~ 1/3 + c cmv A~ 2/3 , with c surf w lOOMeV and c curv w 
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50 100 5 10 15 20 

A A 

Figure 4: As Figure but showing also the liquid drop results. 
300MeV, the stability condition E/A < m n may be written as A > A^ n , where 

^abs = ( C surf + [c^f + 4c curv (m ra - E Q )} l ' 2 V ^ 



Stability at baryon number 30 requires a bulk binding energy in excess of 65 MeV, 
which is barely within reach for m s > lOOMeV if, at the same time, wd-quark matter 
shall be unstable. The proposed cosmic ray strangelet-candidates with baryon num- 



ber 370 |32] would for stability require a bulk binding energy per baryon exceeding 
20 MeV to overcome the combined curvature and surface energies. Absolute stability 
relative to a gas of 56 Fe corresponds to furthermore using 930 MeV instead of m n , 
whereas stability relative to a gas of A-particles (the ultimate limit for formation of 
short-lived strangelets) would correspond to substitution of = 1116MeV. 

Another way of stating the results is to calculate the minimum baryon number 
for which long-lived metastability with respect to neutron emission is possible. This 
requires dE curv /dA + dE sml /dA < m n — e°, or 

^meta _ ( c suri + [ C surf + 3c curv (m n — £ )] \ . . 

min ~ I 3(m n - e°) ) ■ m 



18 



To have A™f* a < 30 requires m n — e° > 30MeV, which is possible, but only for a 
narrow range of parameters. 

This should not, however, defer experimentalists from pursuing the proposed 
searches. After all, the MIT bag model is only an approximation, and in partic- 



ular shell effects can have a stabilizing effect. As stressed by Gilson and Jaffe |22 
the fact that the slope of E/A versus A becomes very steep near magic numbers can 
lead to strangelets that are metastable (stable against single baryon emission) even 
for e° > 930MeV. Also, the time-scale for energetically allowed decays has not been 
calculated. Pauli-blocking is known to delay weak quark conversion in strangelets 



33, 34, Kq, 36|, and this will probably have a significant influence on the lifetimes. 



The existence of small baryon number strangelets is ultimately an experimental issue. 
3.4 Strangelets at finite temperature 

Whereas the calculations above deal with strangelets at zero temperature, the envi- 
ronment in heavy ion collisions is expected to be hot. An advantage of the asymptotic 
mass formula compared to the shell-model calculations is, that it can fairly easily be 
generalized to non-zero temperature. 

The general expression for the thermodynamic potential, Qi, is 

r°° dN- 

a = T9iT / dk— % - In [1 ± exp(-(e(A;) - /i)/T)] (61) 
Jo dk 

where the upper sign is for fermions, the lower for bosons, and the density of states, 
^jjj-, is given by Eq. (|21~1) . For massless quarks (including antiquarks) an integration 
gives, per flavor, 

with a corresponding quark number 



For gluons 



n 9 = -^-T*V + ~T*C. (64) 

The total Q can be found from summing the terms above, and other thermody- 
namical quantities like the free energy and the internal energy can be derived. For 3 
massless quark flavors of equal chemical potential one finds 

11 = { J ^ T ' ~ ^ ~^ + B ) V+ (| T2 + 8^) ° < 65 » 
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E - + V r - + i" 1 + B ) v - + i" 2 ) a < 67 > 

Strangelets are in mechanical equilibrium at fixed temperature and baryon number 
when dF = 0, corresponding to 

In this case one gets the following expressions for the thermodynamic potential, free 
energy, internal energy and baryon number: 

Q = ( il T 2 + C (69) 

V108 4tt 2P J K ' 

.Irpl , 3 4 \ / 41 2 1 o 



P vr 2P / V108 2tt 2P / v ; 

£ = 45V (71) 

Notice that the equations above can also be used in connection with bulk SQM, for 
instance in an astrophysical context, by simply putting C = 0. An external pressure 
can be accommodated by substituting B + P ex temai m place of B. 

Dotted curves in Figure |5| shows the energy per baryon for finite temperature 
strangelets according to the formulae above. Results are given for fixed entropy per 
baryon, where the entropy is calculated from S = — dVt/dT\y^. These results were 
first presented in [RTf. A similar treatment, including finite m s , was published in P§ 



whereas Ref. shows results for a corresponding finite temperature shell model 
calculation, finding that shell structures are washed away at T > 10 MeV, which 
means that liquid drop model and shell model results become indistinguishable at 
high T (S/A). 



As discussed in more detail in [0 further complications arise from the fact, that 



strangelets must be color singlets. This has no influence on the ground state energy 
for T = 0, but for T > quarks are statistically distributed over energy levels, and 
the color singlet constraint reduces the number of possible configurations, forcing the 
energy up for fixed entropy (see also [ID]). The effect is important for A < 100 as 



illustrated in Figure |5|. Similar effects result from insisting that strangelets shall have 
a definite momentum. These destabilising effects can be important in connection 
with experiments, which inevitably create strangelets with rather high entropies. A 
tremendous job remains to be done in calculating the details of strangelet formation, 
evolution, and decay modes, including realistic non-equilibrium effects, etc.! 
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Figure 5: E/A as a function of A for strangelets with equal numbers of massless u, 
d, and s quarks for entropy per baryon of 0, 1, 5 and 10, and B 1 ^ = 145MeV. Solid 
curves include color singlet and zero momentum constraints, dashed curves only the 
color singlet constraint, and dotted curves are without constraints. Entropy increases 
upward. For S — (T = 0) the three curves completely overlap (lowest solid curve). 



4 SQM in cosmology 

4.1 Formation, evaporation and boiling of quark nuggets 

If the cosmological quark-hadron phase transition was first order, supercooling may 
result in concentration of baryon number inside shrinking bubbles of quark phase. 
The amount of baryon concentration depends on the permeability of the "membrane" 
separating the phases and on the turbulent removal of quarks from the phase bound- 
ary. If a quark bubble is able to get rid of entropy fast enough (primarily in the form 
of neutrinos and photons) relative to the rate of baryon number removal, there is a 
chance of reaching baryon number densities in the quark bubbles approaching nu- 
clear matter density. In other words, a quark nugget may form. Whether or not this 
actually happens, or whether one is left with the less extreme, but also interesting 
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scenario where all of the quarks end up in inhomogeneously distributed neutrons and 
protons, giving non-standard Big Bang nucleosynthesis, has been a topic of much 
debate B, HTI, E2|, H3|, E3j , and the final word has probably not been said. 



But even if cosmological quark nuggets do form, they find themselves in a very 
hostile environment with a temperature of order 100 MeV. Under such conditions 
the nuggets are unstable against surface evaporation E5|, [46], [47], H8[ and boiling 



49 , |50|, |yj ; but the crucial question from a cosmological point of view is whether some 



nuggets may survive due to the relatively short time-scale for cooling the Universe 



(the age of the Universe at temperature T being roughly t s 



MeV 



)• 



Alcock and Farhi |45|] showed that the timescale for complete evaporation of a 
quark nugget was smaller than the age of the Universe at temperature T for baryon 
numbers below 

(73) 



4 c~ 

^evap 



2 x 10 56 exp(-3J„/T)/ 1 



3 

n i 



where f n (the phase boundary penetrability of neutrons) was assumed to be close to 
unity. For a homogeneous quark nugget the neutron binding energy /„ = m n — \i u — 2/i^ 
was estimated to be of order 20 MeV. For such a binding, primordial nuggets with 
baryon number A < 10 55 evaporate almost instantly when neutrino heating becomes 
possible at T « 50 MeV pH. 



However, the surface evaporation of neutrons and protons reduces \i u and fid, and 
leads to an increase in fi s . Weak decays, diffusion and convection work to counteract 
this, but the net result is an s-quark enriched layer near the surface. (Small nuggets 
are s-quark enriched throughout their interior). The most efficient way to remove 
the s-quarks is then to emit them in kaons (K°, K~) along with thermal u and d. 
A quasi-equilibrium situation arises with an effective I n ~ 350 MeV [41]]. Thereby 
the baryon number of nuggets surviving evaporation is reduced to 10 46 , and a proper 
inclusion of reabsorption of emitted hadrons (a calculation that has so far not been 
done) may reduce the number somewhat. 

Cosmological nugget evaporation (time-reversed) is closely related to the distilla- 
tion mechanism proposed for strangelet production in relativistic heavy-ion collisions 



57]. There strangeness enhancement occurs due to emission 



3, |53} £0} [LJ, H, |55 
of K+ and K°. 

The calculations described above assume that the penetrability of the phase 
boundary is near 100%. It has been argued that the penetrability may be reduced by 
a few orders of magnitude in a chromoelectric flux tube model. This would decrease 
^evap by a factor permitting smaller nuggets (possibly down to A = 10 39 ) to 
survive 0, [18]]. Again, the limit on A may be further reduced by reabsorption. 

Primordial nuggets are superheated, and may therefore boil by forming bubbles 
of hadronic gas in their interiors [fi5| . However, even though boiling is thermodynam- 
ically allowed, it probably does not play an important role for primordial nuggets (or 
in heavy- ion collisions for that matter), since the time-scale is too short for bubble- 
nucleation to take place |50, . The surface evaporation described above is thus the 



22 



decisive mechanism. 

Some authors have argued ]58|, [59|], that boiling will take place unless a large 
external pressure (e.g. due to a gravitationally bound shell of nucleons) is there to 
prevent it. Such gravitational stabilization only works for masses close to those of 
stars (A ~ 10 57 ). However, the authors discuss only whether boiling is thermo- 
dynamically possible, but neglect that there is not enough time for the bubbles to 
nucleate. 

Apart from trace abundances, one should not expect nuggets smaller than 10 30 — 
10 40 to survive from the early Universe. This however brings one well within the 
causality limit set by the baryon number inside the horizon during the cosmic quark- 
hadron phase transition, 

-.49 



/100MeV\ 2 



A hoi « 10 49 {^ F ^) , (74) 

and includes the "most probable" range of baryon numbers originally predicted by 
Witten ||. It also leaves open the possibility that SQM may explain the dark matter 
problem, and if we understood the details of the quark-hadron phase transition, we 
could even calculate the relative abundances of dark and ordinary matter from first 
principles. 

There is a possibility, that also small traces of primordial nuggets with low baryon 
numbers are left over from the early Universe. Even such traces may in fact be 
"observed" using the astrophysical detectors discussed in Section 0, or via Big Bang 



nucleosynthesis, as explained in Section 4.2 



4.2 Quark nuggets and Big Bang nucleosynthesis 

A crucial property of quark nuggets is the positive electrostatic surface potential of 
the quark phase, which is due to the quarks being stronger bound than the electrons 
(electrostatic forces are weaker than strong forces). For typical nugget parameters 
the electrostatic potential can be several MeV, so except at very high temperatures, 
protons and nuclei are repelled from nuggets, whereas neutrons are absorbed, adding 
one unit of baryon number. 



This opens the intriguing possibility of using SQM as an energy source at 
least in principle. It also makes it possible to use Big Bang nucleosynthesis as well 
as the properties of pulsars to place very stringent limits on the abundance of quark 
nuggets in the Universe. 

During Big Bang nucleosynthesis (T < 1 MeV), nuggets absorb neutrons but 
not protons. This means that the presence of quark nuggets reduces the neutron- 
to-proton ratio, thereby lowering the production of 4 He. The helium-production 
is very sensitive to the total amount of nugget-area present, and in order not to 



ruin the concordance with observations, one finds [|60| that only nuggets with A > 



Abbn ~ 10 23 fi n /i 6 / 3 are allowed during nucleosynthesis. Here f2 nug is the present- 
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day nugget contribution to the cosmic density (in units of the critical density) , h is the 
Hubble parameter in units of 100 km sec -1 Mpc -1 , and f n < 1 is the penetrability 
of the nugget surface. Slightly stronger limits can be obtained from inclusion of 
all light nuclei instead of 4 He only |6lJ . (Ref . |62| found good correspondence with 
nucleosynthesis for a nugget-dominated, Q — 1 Universe if A « 10 17 , but as shown in 
[ [46[1 , this was due to an erroneous emission rate for nucleons.) 

The nucleosynthesis calculations leading to A BBN neglected inhomogeneities in the 
nucleon distribution, and all nuggets were assumed to have the same baryon number. 
However, the formation of 4 He is an on-off process over a limited range of A, so the 
detailed behavior of the inhomogeneities may not be so important. 

Note that SQM, in spite of it carrying baryon number, does not contribute to the 
usual nucleosynthesis limit on i\ aryon . The SQM baryon number is "hidden" in quark 
nuggets long before Big Bang nucleosynthesis begins, and the nuggets only influence 
nucleosynthesis if they have a big total surface described above. 

Evaporating nuggets would lead to strongly inhomogeneous nucleosynthesis with 
enhanced heavy-element formation. This aspect has recently been studied in p3l . 



4.3 Quark nuggets as dark matter 

Witten argued that quark nuggets might be a natural explanation of the cosmo- 
logical dark matter problem, in principle allowing a calculation of the relative amount 
of dark matter and ordinary baryons. In view of the evaporation discussed above, 
this idea now seems less likely, but is certainly not ruled out for A > 10 30 . Massive 
quark nuggets decouple from thermal equilibrium with the radiation bath very early 
in the history of the Universe, quickly slow down, and behave as cold dark matter in 
the context of galaxy formation. 

Of course it should again be noted, that all of the interesting cosmological conse- 
quences of the quark-hadron phase transition require the transition to be first order, 
in agreement with recent lattice QCD calculations. 



5 SQM in neutron stars; strange stars 

It has been known for many years, that neutron stars may in fact be "hybrid stars" 
consisting of "ordinary" nuclear matter in the outer parts and quark matter in the 
central regions. This will be the case if SQM is metastable at zero pressure, being 
stabilized relative to hadronic matter by the high pressure within a neutron star 

II Hill- 

If SQM is absolutely stable at zero pressure, an even more intriguing possibility 



opens up, namely the existence of "strange stars" g, |67|, gy] consisting completely 
of SQM (perhaps apart from a minor crust to be discussed below). Such strange stars 
behave quite differently from neutron stars due to the unusual equation of state. For 
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massless quarks the total energy density is given by p = p q + B, and the total pressure 
by P = P q — B, where p q is the energy density of quarks, and the pressure of the 
quarks is P q = p q /3, since massless quarks are relativistic. The equation of state is 
thus given by 

P= l -{p-AB). (75) 
The exact equation of state taking into account m s ^ is very similar [68] since 



s-quarks are relativistic for low m s and not present for high m s . (I here assume that 



a s = 0, but recall from Section |2.3| , that a non-zero a s effectively corresponds to a 
reduction of B). 

The structure of a strange star is calculated from the Oppenheimer-Volkoff equa- 
tion, describing the balance between gravity and pressure gradient, using the equation 
of state given above. The surface of the star corresponds to P = 0, a condition ful- 
filled for p = 4B, which for typical values of B is somewhat more than the density of 
ordinary nuclear matter! For stellar masses below 1M (M is the solar mass) this 
density is almost constant throughout the star, so to a good approximation total mass 
and radius are related by M cx R 3 , a relation in striking contrast to ordinary neutron 
stars, where M cx R~ 3 . This means that low-mass neutron stars and strange stars 
have widely different radii, possibly allowing observational distinction. Unfortunately 
Nature prefers to form these compact objects with masses near 1.4M Q , according to 
stellar evolution models. For such a mass gravity rather than bag pressure plays the 
dominant stabilizing role, and there is no significant difference between neutron star 
and strange star radii. Also the maximum mass given by gravitational instability 
(the Chandrasekhar limit) is similar, of order 2M Q . In contrast to ordinary neutron 
stars, which are unstable for masses below 0.1 M Q , strange stars have no minimum 
mass; the sequence continues smoothly to the domain of strangelets. 

For the simple equation of state discussed above, the only natural energy scale in 
the problem is B 1 ^. Thus there exists a homology transformation between strange 
star models for different values of B. In particular, the maximum mass of a strange 
star is given by 

M max = 2.006B U \ /2 M Q . (76) 

The corresponding minimal radius, maximal moment of inertia, maximal central 
density, surface density, and minimal rotation period (the so-called Kepler period 
corresponding to mass-shedding at the equator), are given by 

R min = lO.Q^fkm, (77) 

J max = 2.256 x lO^fg cm 2 , (78) 
/W = 1.97 x 10 1B Bi4sg cm" 3 , (79) 
p surf = 4.102 x 10 14 5 145 g cm" 3 , (80) 
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P min = OMBi£( 2 ms. (81) 

Bare strange stars (strange stars with quark matter all the way to the surface) 
have quite unusual properties. The density abruptly jumps from to p sur f (Eq. (|80|) ), 
and the density is almost constant through the interior (except when the mass is 
close to M max ). The plasma frequency of the star is huge, meaning that photons 
with energies below 20MeV are reflected from the surface, whereas the star itself 



can only emit photons with higher energies [68, 70). Even more important, because 



of the strong interaction binding of the surface material, the star is not subject to 
the "Eddington limit", which for ordinary neutron stars limits the luminosity to be 
below 10 38 erg/s (for higher luminosities the radiation pressure would exceed the 
gravitational attraction and expel the surface layers). As discussed below, this could 
lead to important "applications" of strange stars. 

This approach may however be oversimplified because real strange stars may have 
surfaces more like ordinary neutron stars. In particular, a solid crust of ordinary 
material may form from accretion by the strange star after formation, or from material 
that was not converted during neutron star burning (see Section |5~5"1) . Such a crust 
may be held up by the extreme, outward directed electrostatic potential of 10 17 - 
10 18 V/cm, created by the electron atmosphere with a thickness of a few hundred 
Fermi. This atmosphere merely expresses that the electrostatic binding of electrons 
is weaker than the strong binding of quarks; therefore the electron distribution does 
not end abruptly like that of quarks (the detailed structure was found from a Thomas- 
Fermi calculation by Alcock, Farhi and Olinto see also ||71[|). 

The electrostatic potential can sustain a significant crust of ordinary neutron star 
material. The limit is given by the neutron drip density (4 x 10 11 gem -3 ), above 
which neutrons drip out of nuclei and would be swallowed by the quark phase. This 
crust may be decisive for interpretation of pulsar behavior (Section |5.2p . 

As emphasized by Glendenning, Kettner and Weber [72|, [73[] , the existence of 
crusts not only changes the mass-radius relation for strange stars, but also opens a 
rich plethora of new stellar configurations. In particular, one may have a sequence of 
"strange dwarfs", much like white dwarfs except for an SQM core. At present there 
is no well-studied model for formation of such strange dwarfs. 

Another possibility for formation of a (solid?) crust has been suggested in f74| . 
This mechanism relies on the existence of stable, low-baryon number strangelets (in 
this context sometimes denoted "quark-alphas" for the A = 6 strangelet analog of a 
helium nucleus [[7511 ) which could act as "nuclei" in the surface region. Whereas this 
possibility may seem less likely from the discussion in Section |3|, it can not be entirely 
ruled out. 

Finally, it is worth noticing that Glendenning [[T!J has argued that neutron stars 
may contain regions with mixed quark and hadron phases. (This possibility was 
missed in earlier studies due to an erroneous assumption of local rather than global 
charge neutrality). Depending on parameters the mixed phase region can occupy 
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a significant fraction of the star, and may show unusual topologies (plate-like or 
cylinder-like structures, rather than just spherical quark bubbles embedded in hadrons 
or vice versa W% [F3, 



Studies of strange stars have not been pursued to the degree of detail known for 
ordinary neutron stars, and it is premature to draw any detailed conclusions. How- 
ever, in the following, I shall look at some of the properties expected and emphasize 
the possible observable differences between neutron stars and strange stars. 

5.1 Neutrino cooling 

A distinction between strange stars and neutron stars was for a long time believed 
to be a much more rapid cooling of SQM due to neutrino emitting weak interactions 
involving the quarks f68|. Thus a strange star was presumed to be much colder than 



a neutron star of similar age, a signature potentially observable from x-ray satel- 
lites. Only a few speculative mechanisms, such as the existence of kaon condensates 
might mimic the speed of quark matter neutrino cooling. Recently the story has 
been complicated considerably by the finding that ordinary neutron /9-decay may be 
energetically allowed in nuclear matter [80|, so that the cooling rate can be compa- 



rable to that of SQM. For this reason I shall not discuss the issue here, but refer the 
reader to an excellent review of neutron star cooling by Pethick |pl |, and a recent 
reinvestigation of strange star cooling by Schaab et al. |82| . 



5.2 Pulsar glitches 

One important feature seems to distinguish strange stars from neutron stars in a 
manner with observable consequences, and that is the distribution of the moment of 
inertia inside the star. Ordinary neutron stars older than a few months have a crust 
made of a crystal lattice or an ordered inhomogeneous medium reaching from the 
surface down to regions with density 2 x 10 14 gcm -3 . This crust contains about 1% 
of the total moment of inertia. Strange stars in contrast can only support a crust 
with density below the neutron drip density (4.3 x 10 11 gcm~ 3 ). This is because free 
neutrons would be absorbed and converted by the strange matter. Such a strange 
star crust contains at most a few times 10~ 5 of the total moment of inertia. This is 
an upper bound, since the strange star may have no crust at all, depending on its 
prior evolution. And recent studies of the mechanical balance between electric and 
gravitational forces on the crust indicate, that only densities up to perhaps 10 n g cm -3 
may be achieved J83], j84| . 

As stressed by Alpar ||851 , and also pointed out by Haensel, Zdunik, and Scha- 
effer JF/J, and by Alcock, Farhi, and Olinto ||68|| , this difference in the moment of 



inertia stored in the crust of neutron stars and strange stars seems to pose sig- 
nificant difficulties for explaining the glitch-phenomenon observed in radio pulsars 
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with models based on strange stars. Glitches are observed as a sudden speed- 
up in the rotation rate of pulsars. The fractional change in rotation rate Q is 
AQ/Q w 1CT 6 — 1CT 9 , and the corresponding fractional change in the spin-down 
rate Q is of order Ati/Cl w 10~ 2 — 10~ 3 . Regardless of the detailed model for the 
glitch phenomenon these jumps must involve the decoupling and recoupling of a 
component in the star containing a significant fraction, /$//, of the total moment 
of inertia; fk/I = AI/I w Att/Q « 10^ 6 — 10~ 9 (Alpar actually argued that 
fli/I ~ Ati/Cl pa 10~ 2 — 10~ 3 , where / is the fractional change in J^, but this is not 
necessary ||36| ). This role is played by the inner crust of an ordinary neutron star, but 
the crust around a strange star is smaller; less than a few times 1O _5 M with I CIUS t/I 
around a few times 10 -5 for ordinary neutron star masses of 1.4M (higher for less 



massive stars). These numbers are based on models by Glendenning and Weber [86 
assuming a maximum mass crust, i.e. a crust reaching neutron drip density at the 
base, so it seems fair to conclude, that strange stars in fact may have sufficiently 
massive crusts to account for glitches, but that parameters in that case are fairly 
tightly constrained. 

Other possibilities for glitches in strange stars could involve a crust composed of 
strangelets (cf. the "quark-alpha" scenario in @), not to mention the possibility of 



a quark-hadron mixed phase flTBL 77, 78, 79]. There is still a lack of any detailed 



model for how the magnetic field structure and other crucial aspects of a pulsar 
can be modeled for strange stars. Presumably a strange star cannot do the job 
without significant structure, such as a crust and/or superfluidity/superconductivity 
in certain regions. These issues have only been very superficially studied and need 
further consideration. The present lack of such models should not be used to dismiss 
the possibility of strange stars. 

5.3 Strange star oscillation and maximum rotation rate 

One of the most interesting differences between neutron stars and strange stars is 
related to the damping of instabilities. 

First it should be noticed that a strange star is a very stable system. Strange 
stars may have radial oscillations with a fundamental period of 0.06-0.3 ms [j57||, but 



these are characterized by rapid damping in a matter of seconds |8q , p9j . p0| , |9T| . This 
is due to the extremely high viscosity of SQM. 

The large viscosity also plays a role in setting the maximum rotation limit for 
strange pulsars (or hybrid stars with SQM cores). The ultimate rotation limit corre- 
sponds to mass-shedding from the equator of the star (this is called the Kepler limit 



and is of order 0.6 msec for a strange star, Eq. (plf) ; see Zdunik |92j for a review). 
But before reaching such rotation rates, the pulsars become unstable to non-radial 
deformations and are slowed down by emission of gravitational radiation. Shear and 



bulk viscosities tend to stabilize the star against these instabilities [SDL 031], and the 
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high value for the bulk viscosity may mean that strange pulsars in contrast to ordi- 
nary pulsars can reach submillisecond periods ||91|| . Thus the discovery of very fast 



pulsars may be an indication favoring the existence of strange stars. 

And even more exciting, it has been shown over the last few months, that ordinary 
neutron stars when they are young and hot are subject to a new class of instabilities, 
called r-mode instabilities |94|, [5^, j96| , which during their first year of existence slows 



the rotation rate to only a few per cent of the Kepler limit. Rotation periods faster 
than 10 msec are unlikely after that, until some pulsars at a much later age may be 
spun-up by angular momentum transfer in binary systems, and thereby explain the 
rapid old pulsars with periods down to 1.56 msec. In contrast, strange stars are not 
subject to these instabilities until they are thousands of years old, and even then 
only for periods faster than 2-3 msec ||97|| . This seems to imply, that the most robust 
signature for the existence of strange stars (or neutron stars with a substantial fraction 
of high viscosity quark matter in the interior) is to search for young pulsars with 
rotation periods below, say, 5 msec (even stars with longer periods may candidate). 
These can not be ordinary neutron stars, whereas quark matter is the only substance 
known to have a bulk viscosity high enough to offer an explanation. 

The bulk viscosity of strange quark matter depends on the rate of the non-leptonic 
interaction 

u + d «-> s + u. (82) 



(The rate for this reaction has recently been calculated by Madsen |35 , and Heiselberg 
[ |36|1 ; earlier studies, including that of Ref. are incorrect). This reaction changes 
the concentrations of down and strange quarks in response to the density changes 
involved in vibration or rotational instabilities, thereby causing dissipation. This 
dissipation is most efficient if the rate of reaction fl32|) is comparable to the frequency 
of the density change. If the weak rate is very small, the quark concentrations keep 
their original values in spite of a periodic density fluctuation, whereas a very high 
weak rate means that the matter immediately adjusts to follow the true equilibrium 
values reversibly. But in the intermediate range dissipation due to PdV-work is 
important. 

The importance of dissipation due to Eq. ( |8"2"D was first stressed by Wang and Lu 
[|S8|j in the case of neutron stars with quark cores. These authors made a numerical 
study of the evolution of the vibrational energy of a neutron star with an O.2M quark 
core, governed by the energy dissipation due to Eq. (|82|). Sawyer |89|] expressed the 
damping in terms of the bulk viscosity, a function of temperature and oscillation 
frequency, which he tabulated for a range of densities and strange quark masses. 
Sawyer's tabulation has later been used in studies of quark star vibration and 
of the gravitational radiation reaction instability determining the maximum rotation 
rate of pulsars [p3| . The latter study concluded, that the bulk viscosity is large 



enough to be important for temperatures exceeding 0.01 MeV, but that it should be 
a few orders of magnitude larger to generally dominate the stability properties. 
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However, as has been pointed out in [pi]] , the bulk viscosities calculated in |89 
depend on the assumption, that the rate of Eq. (|82f) can be expanded to first order 
in djj, — jMs — Pdi where /ij ~ 300MeV are the quark chemical potentials. This 
assumption is not correct at low temperatures (2irT -C 5fi), where the dominating 
term in the rate is proportional to 5/i 3 . Furthermore, the rate in [^] is too small 
by an overall factor of 3, and a discrepancy of 2-3 orders of magnitude, perhaps 
due to unit conversions, appears as well. Taken together, these effects lead to an 
upward correction of the bulk viscosity by several orders of magnitude, and thereby 
increases the importance for the astrophysical applications. The non-linearity of the 
rate also means, that the bulk viscosity is no longer independent of the amplitude of 
the density variations. The resulting bulk viscosity is (in cgs-units, with m s , T, and 



fi d « 235MeV(p/p B 



a/3 



C « 3.09 x 10 28 m> 



in MeV, and the oscillation frequency w in s *) 



3 / m 2 Av * 



4 V 3fi d v 



+ 4vr 2 T 2 



gem 1 s . 



(83) 



For typical values 



100 MeV, 



/V 



300MeV, 



2 x 10 4 s- 1 ) this is C 



1.6 x 10 28 [93(Av/v ) 2 + 39T 2 ] gcm _1 s , where Av/vq is the perturbation amplitude. 

For a star of constant density (an excellent approximation for a strange star, 
except very close to the gravitational instability limit) Sawyer |[{J estimated the 
damping time as 

t d ss 1.5 x 10 25 C _1 s. (84) 

Thus, even at very low temperatures, high amplitude oscillations are damped in 
fractions of a second, and those of low amplitude in a matter of minutes, if one takes 
into account, that the temperature of the star increases due to the heat released by 
viscous dissipation, which can speed up the damping of vibrations. 



The discussion above was based on rather crude estimates |91|| . A detailed, general 
relativistic, numerical treatment along the lines of Cutler et al. |K| is clearly needed. 

As mentioned previously, viscosity also plays an important role in setting the 
maximum rotation rate of pulsars. Gravitational radiation reaction instabilities (as 
opposed to "Keplerian mass- shedding" ) is supposed to set the ultimate rotation rate 
limit, but the larger the damping by shear and bulk viscosity is, the closer the rate 
can get to the Keplerian limit given in Eq. (|8~If) . 

The shear viscosity of SQM due to quark scattering has recently been recalculated 
by Heiselberg and Pethick fl98[. Their results for T /x can be written as 



, 1\ 5 / 3 

15 / V- 1 \ 



r] w 4.0 x 10 iO — 




gem *s . 



Investigations by Colpi and Miller |K| based on the older viscosities in 
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55) 



in- 



dicated, that the minimal rotation period of strange stars might be set by the grav- 
itational radiation reaction instability of m = 2 or m = 3 modes at or just below 
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1 millisecond. With the new, much larger, viscosities, the non-axisymmetric insta- 
bilities will be suppressed, and it is not unreasonable to expect, that the maximum 
rotation frequency of strange stars will be close to the Keplerian limit. Detailed 
numerical calculations like those in Colpi and Miller [f?3j| , including the new viscosi- 
ties and effects of dissipative heating, are required to settle the issue, but they are 
complicated by the non-linear behavior of the new bulk viscosity. 

Whether or not the ultimate rotation period of strange stars can be significantly 
smaller than for neutron stars is of importance for old pulsars spun-up by accretion. 
But perhaps the most clear-cut signature for the existence of strange stars would be 
the (almost) lack of sensitivity to r-mode instabilities, which as mentioned earlier 



allows young strange stars to rotate much faster than young neutron stars [97 



5.4 Gamma-ray bursters 

Strange stars because of their high surface density, strong binding (making it pos- 
sible to circumvent the Eddington limit), and special emission properties have been 
suggested as explanations for some of the more mysterious cosmic events, namely 
7-ray bursters. These are bursts of 7-rays of a few seconds duration, coming from 
unidentified sources which are presumably at extragalactic distances. 

No consensus exists concerning the nature of these bursts, but Alcock, Farhi and 



Olinto ||100|| suggested a detailed model for the most prominent of the bursters, the 
one on 5 March 1979. Their model is based on an impact of a 1O _8 M lump of 
SQM on a rotating strange star, and the authors are able to explain most of the 
observations concerning energetics and time-scales under the assumption that the 
burster is located in a supernova remnant in the Large Magellanic Cloud, as position 
measurements seem to indicate. An alternative model for this source and for soft 
7-repeaters in the framework of strange stars with "quark-alpha" surface properties 



was suggested in [ 101 1 . Other strange star models for soft 7-repeaters and x-ray 



bursters include [|102|, |103f 



7-ray bursters at truly cosmological distances could be due to collisions of two 
strange stars in binary systems ||104j| , each collision releasing 10 50 ergs in the form of 



gamma rays over a time-scale of 0.2 s. 

There are, however, literally hundreds of different models for 7-bursts, and in 
spite of improved observational data the interpretation is at present unclear. 

A recent identification of the x-ray source Her X-l as a strange star [|105|| was 
unfortunately based on incorrect use of bag model parameters [ |106|| . 

5.5 Formation of strange stars 

If strange quark matter is stable, strange stars may be formed during supernova- 
explosions, and neutron stars can be converted to strange stars by a number of 
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different mechanisms, such as pressure-induced transformation to uds-quark matter 
via ud-quark matter, sparking by high-energy neutrinos, or triggering due to the 
intrusion of a quark nugget. These and other possibilities were described by Alcock, 
Farhi, and Olinto . 

As soon as a lump of strange matter comes in contact with free neutrons it starts 
converting them into strange matter. The burning of a neutron star into a strange 



star was discussed by Baym et al. [ 107|| and Olinto ||108| , and it was shown that 



the star would be converted on a rather small time-scale set by quark diffusion and 
flavor-changing weak interactions. (The huge difference in the speed of the conversion 
front found in these papers is partly due to the omission of a factor c 1 / 2 , where c is the 
speed of light, in equation (6) of Baym et al.) Later studies ||109| , |110|| found burning 



times in the range of 1-10 seconds under various parameter assumptions (see also 



Olinto [l I l] for a review). For the fastest burning times, the energy liberated may 



be important for the supernova mechanism and supernova neutrino bursts. Horvath 



and Benvenuto ||112|| have questioned the stability of "slow" neutron combustion 
and suggested that the conversion takes place much faster as a detonation. So far, 
the investigations of neutron star burning have been rather crude, neglecting many 
aspects of transport theory, heat conduction etc. A detailed study of this phenomenon 
would be interesting. 

Perhaps the most likely mechanisms for initiating the formation of a strange star 
involves either a seed of SQM in the star (see Section ^|), or thermal formation of 
quark matter bubbles. Thermal triggering of neutron star transformation may be 
understood qualitatively in terms of simple boiling theory. Before considering a more 
realistic equation of state it is instructive to study the boiling of a pure neutron 
gas into quarks. The quark bubbles formed consist of u and d quarks in the ratio 
1:2; only later weak interactions may change the composition to an energetically 
more favorable state. Thus quark chemical potentials are related by fid = 2 l ^fjb u , 
and fi n = fi u + 2 = (1 + 2 4 / 3 )/i M , assuming chemical equilibrium across the phase 
boundary. 

The free energy involved in formation of a spherical quark bubble of radius R and 
volume V is given by 

F = -APV + 8tt 7j R (86) 

where 

AP = P ud -P n = 4^ ~B-P n (87) 
and the curvature energy coefficient 

7 = -g^- (88) 
The free energy has a maximum at the critical radius 

r c = (27/AP) 1 / 2 (89) 
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and the corresponding free energy 



W c = F{r c ) = 167r 7 r c /3 (90) 

is the work required to form a bubble of this radius which is the smallest bubble 
capable of growing. It is a standard assumption in the theory of bubble nucleation 
in first order phase transitions that bubbles form at this particular radius at a rate 
given byQ 

K « T 4 exp(-W c /T). (91) 

The simplest possible equation of state for the neutron gas is that of a zero 
temperature, nonrelativistic degenerate Fermi-gas, where 

Pn = ^" T (92) 
lo7T m„ 



and the baryon density 

3tt 2 

A necessary condition for boiling is that AP > 0. This leads to an upper limit 
on the bag constant, -B max , from Eq. (|87|) as illustrated in Fig. || (from ||113|| — the 



n B = ^ - (93) 



corresponding limit for the Bethe- Johnson equation of state is shown for comparison; 
it is seen to be very similar). This was also used as a criterion for neutron star 
stability by Krivoruchenko and Martemyanov ||11 1 | . 



Also shown in Fig. |6| is the limit on the bag constant below which bubble nucle- 
ation takes place at rates exceeding 1 km _3 Gyr _1 and 1 cm~ 3 s _1 , respectively, for 
temperatures of 1, 2, 3 and 10 MeV (-B m ax can be considered as the limit for infinite 
temperature). One notes that the possibility of bubble nucleation is fairly insensitive 
to the temperature as soon as T exceeds a few MeV, whereas thermally induced bub- 
ble nucleation is impossible for T < 2MeV (recall from Section |2.2| that the stability 
of ordinary nuclei against decay into quark matter requires that B > (145MeV) 4 ). 
This confirms an estimate in [|115|| (see also [ 116 , |1 1 7|| ) . The range of bag con- 



stants for which a hot neutron star may transform into quark matter is thus roughly 
145MeV < B 1 ^ < 155MeV. 

Results for a more realistic mean field equation of state are also shown in the 



Figure (see ||113|1 for further details). While the detailed numbers change, the over- 
all conclusion does not. Quark matter bubbles may nucleate (possibly followed by 
burning of the star into SQM) in neutron stars/supernovae if the bag constant is low, 
and if the temperature exceeds a few MeV (thus the process is most likely during 
the supernova explosion itself). Should thermal nucleation not take place, one of 
the other mechanisms mentioned above must be relied on. Apart from seed-induced 
burning, all of these are likely to be much less efficient than thermal nucleation. 

2 The prefactor may differ from T 4 , but this is of minor practical importance due to the dominant 
exponential. 
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Figure 6: Upper limits on the bag constant allowing thermal nucleation of quark 
matter bubbles in neutron stars as a function of baryon number density in the hadron 
phase. To the left are shown results for the simple neutron gas model discussed in 
the text. To the right for a more realistic mean field approximation. See text for 
further explanations. 



6 SQM in cosmic rays 



De Rujula and Glashow [ ITS] argued that unusual meteor-events, earth-quakes, etched 
tracks in old mica, in meteorites and in cosmic-ray detectors might be used for ob- 
servation of quark nuggets hitting the Earth or its atmosphere. In particular they 
were interested in the possibility of detecting a galactic dark matter halo of nuggets, 
where typical velocities would be a few hundred kilometers per second, given by the 
depth of the gravitational potential. Even if nuggets only survived from the Big 
Bang in small numbers, or were spread in our galaxy by secondary processes such as 
strange star collisions, there could be a potentially observable flux of nuggets hitting 
the Earth. The only data actually investigated in their paper came from a negative 
search for tracks in ancient mica, and corresponded to a lower nugget flux limit of 
8 x 10 -19 cm -2 s^ 1 sr _1 , for nuggets with A > 1.4 x 10 14 (smaller nuggets would be 
trapped in layers above the mica samples studied). This can also be expressed as an 
excluded range of 1.4 x 10 14 < A < 8 x 10 23 p 2 4^250; where v = 250km s~ 1 i7 25 o and 
p = lCr 24 gcm~ 3 p 2 4 are the typical speeds and mass density of nuggets in the galactic 
halo. The speed is given by the depth of the gravitational potential of our galaxy, 
whereas P24 ~ 1 corresponds to the density of dark matter. In these units the number 
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of nuggets hitting the Earth per cm 2 per second per steradian is 6.0 x 10 5 A~ 1 p 2 4f25o- 
Later investigations have improved these flux limits somewhat. These Earth- 
based flux-limits ||119| , 12Q| are shown in Figure [?|. It is seen that quark nuggets with 



3 x 10 7 < A < 5 x 10 25 seem incapable of explaining the dark halo around our galaxy, 
but a low flux either left over from the Big Bang or arising from collision of strange 
stars cannot be ruled out. If the strange matter hypothesis is valid, one should indeed 
expect a significant background flux from stellar collisions, since several pulsars are 
members of binary systems, where the two components are ultimately going to collide. 
If such collisions spread as little as O.1M of non-relativistic strangelets with baryon 
number A, a single collision will lead to a flux of 10 _6 A _1 f25o cm~ 2 s~ 1 sterad _1 , 
assuming strangelets to be spread homogeneously in a halo of radius 10 kpc. 



Such a flux-level is below the sensitivity of present experiments, but Madsen ||121 
suggested that neutron stars and their stellar "parents" may be used as alternative 
large surface area, long integration time detectors. The reason is simple. The presence 
of a single quark nugget in the interior of a neutron star is sufficient to initiate a 
transformation of the star into a strange star [EL KSL |107||. The time-scale for the 



transformation is short, between seconds and minutes ||107| |108| , |111| , |109| , |110]| , so 
observed pulsars would have been converted long ago, if their stellar progenitors 
ever captured a quark nugget, or if the neutron stars themselves absorbed one after 
formation. 

The rate at which quark nuggets hit the surface of a star depends on the phase 
space distribution of nuggets relative to the star. For an infinite bath of positive 
energy nuggets with an isotropic, monoenergetic distribution function, the number 
accretion rate is given by 

F = 1.39 x lO 3 ^- 1 A' 1 MRp 24 V25o [l + O.WAv^RM- 1 ] , (94) 

where M and R denote the stellar mass and radius in solar units. For the Sun 
the second term in parenthesis (the geometrical term) contributes only slightly to the 
accretion rate, and the contribution is even less important for more massive stars and 
for compact objects like white dwarfs and neutron stars (in contrast, the geometrical 
term dominates for accretion onto the Earth). In the following I therefore only take 
the first term (gravitational) into account. 

To convert a neutron star into strange matter a quark nugget should not only 
hit a supernova progenitor but also be caught in the core. Similarly, nuggets hitting 
a neutron star after its creation have to penetrate the outer layers and reach the 
neutron drip region. These issues were discussed in ||121]| . 

A main sequence star is capable of capturing quark nuggets with baryon numbers 
below Astop, where 

A STOP = 5.0 x 10 31 M~ L8 . (95) 

This works for non-relativistic nuggets, which are basically braked by inertia, i. e. 
they are slowed down by electrostatic scatterings after plowing through a column of 
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mass similar to their own, and afterwards settle in the stellar core. In particular it is 
valid for nuggets moving with virial speed in our galactic halo. Relativistic nuggets, 
like those reported in some cosmic ray observations, may be destroyed after collisions 
with nuclei in the stellar atmospheres, and so the limit can not be used immediately, 
but it is worth noticing, that even a tiny fraction of a nugget surviving such an event 
and settling in the star is sufficient to convert the neutron star to a strange star. 
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Figure 7: Astrophysical flux-limits [|121|| compared to the flux expected for a galactic 
halo of nuggets being the dark matter [ 118|| , and to the experimental results for cosmic 
rays ||1 19| , |120|| . The three horizontal parts of the solid curve correspond to capture 



in main sequence supernova progenitors, post main sequence stars, and neutron stars 
younger than the Vela pulsar (10 4 years). 

For nuggets with A < A^top the sensitivity of main sequence stars as detectors 
is remarkable, as it is given by the limit of one nugget hitting the surface of the 
supernova progenitor in its main sequence lifetime! Converted into a flux, J 7 , of 
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nuggets hitting the Earth per cm 2 per sec per steradian, it corresponds to 



f = 4x l(T 42 M ai *4 . (96) 

As can be seen from Figure [7], this is a factor of 10 20 -10 40 more sensitive than ordinary 
experiments! 

If it is possible to prove that some neutron stars are indeed neutron stars rather 
than strange stars, the sensitivity of the astrophysical detectors rules out quark 
nuggets as being the dark matter for baryon numbers in the range A < 10 34-38 . 
And it questions the whole idea of stable strange quark matter, since it seems im- 
possible to avoid polluting the interstellar medium with nuggets from strange star 
collisions or supernova explosions at fluxes many orders of magnitude above the limit 
measurable in this way. 

If on the other hand SQM is stable, then all neutron stars are likely to be strange 
stars, again because some pollution can not be avoided. 

The Sun would in this way accrete 3.7 x iO~ 2O p 2 4'f ; 2~5o^0/y ear > or a total of 
iO~ lo p24t>25o^0 in its total lifetime on the main sequence. Very low-mass nuggets 
collected near the solar center in this manner might have an impact on the energy 
production ||122|| , but the effect is negligible unless the electrostatic barrier at the 



nugget surface is much smaller than expected, or unless very special circumstances 
allow nuggets to catalyze nuclear reactions |21] . 



The Sun will develop into a white dwarf in about 6 x 10 9 years. As just mentioned, 
the Sun would accrete a core of lO -10 /^^^^© m its total lifetime on the main se- 
quence. Such accretion is too small to lead to a strange dwarf distinguishable from 
an ordinary white dwarf as suggested in [72 , 73|. However, higher concentrations 



could occur if quark nuggets were somehow mixed into the gas cloud from which the 
star originally formed. Whether this is likely to happen depends strongly on assump- 
tions regarding the velocity distribution of the nuggets formed, and the possibility of 
interactions with the gas ||123| . 



Most of the discussion above dealt with halo nuggets moving at non-relativistic 
velocities. Relativistic nuggets are not as easily detected using neutron stars, since 
they may be destroyed in collisions with nuclei in the star. On the other hand two 
relativistic candidate events with charge Z = 14 and mass A « 370 were found in a 
balloon experiment by Saito et al. [32]. This corresponds to a rather high flux, and 



it is not quite clear how to produce such nuggets, though spallation of larger nuggets 
originating from strange star collisions may be involved [|1 24|] . 

Quark nuggets have also been suggested as candidates for the Centauro cosmic-ray 
events |2|, |125| , |126|| . Centauro primaries may have a flux as high as 10~ 14 cm~ 2 s _1 and 



A ~ 10 3 . Since Centauro primaries move at relativistic speeds they are destroyed by 
inelastic collisions when hitting a star, so the flux-limits given above cannot directly 
be used to rule out quark nuggets as Centauro primaries. However the mechanism 
producing the primaries must be tuned so that it only produces relativistic quark 
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nuggets in order not to conflict with the flux-limits in Figure |^ for non-relativistic 
nuggets. 



7 Conclusion and outlook 

The possible stability of strange quark matter is a fundamentally exciting idea. 
Should it turn out to be true, many textbooks in nuclear, particle and astrophysics 
will need revisions, but our daily lives will not be changed dramatically, apart from 
possible technological applications such as energy production and disposal of radioac- 



tive waste [14], |127 . There are two main reasons why stability of SQM is possible 
without drastic consequences. The first reason is that stability requires a certain 
minimum strangeness content, so ordinary nuclei do not decay into strangelets. The 
second reason is the positive electrostatic potential of the quark phase in a strangelet, 
which means that you could walk around with a lump of SQM in your pocket without 
being swallowed. 

While heavy- ion collisions is the way to look for small (meta) stable strangelets, 
astrophysics gives a possibility for testing larger (and therefore more stable) SQM- 
systems. Direct cosmic ray searches is an obvious way to look for intermediate baryon 
numbers in the form of relativistic or non-relativistic lumps produced in strange star 
collisions, and for leftovers from the Big Bang. The latter can only exist for very high 
baryon numbers (cf. Section whereas a galactic background of the former seems 
unavoidable if the strange matter hypothesis is correct. 

Strange stars may be the most promising place to look for SQM, but as explained 
in Section |^ it is actually hard to find clear-cut ways of distinguishing strange stars 
from neutron stars, unless one finds an object of very low mass. Pulsar rotation 
properties at present seem to provide the best clue, in particular after the finding 
that young strange stars in contrast to neutron stars are not braked by gravitational 
wave emission due to r-mode instabilities. 

If SQM is only metastable, heavy-ion physicists may still have a chance of finding 
it; the cosmo logical quark-hadron phase transition may still lead to inhomogeneities 
of importance for Big Bang nucleosynthesis (without quark nuggets left over); and 
neutron stars may still have strange matter cores. 

In any case the confirmation or disproof of the existence of (meta)stable strange 
quark matter via experiments and astrophysical tests makes it possible to limit strong 
interaction parameters that are otherwise difficult to probe. This in itself is a good 
reason for continued studies of the physics and astrophysics of strange quark matter. 
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